Test of finite temperature RPA on a Lipkin model 
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We investigate the applicability of finite temperature random phase approximation (RPA) using 
a solvable Lipkin model. We show that the finite temperature RPA reproduces reasonably well the 
temperature dependence of total strength, both for the positive energy (i.e., the excitation) and the 
negative energy (i.e., the de-excitation) parts. This is the case even at very low temperatures, which 
may be relevant to astrophysical purposes. 
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The random phase approximation (RPA) [1] and its 
extension, quasi-particle RPA (QRPA), have successfully 
described nuclear giant resonances[2, 3], Usually, those 
giant resonances are built on the ground state, but the 
collective states can be excited also from excited states 
[4]. In the early 80's, heavy-ion fusion experiments re- 
vealed the existence of giant resonances in hot nuclei [5— 
9], and properties of collective excitation at finite tem- 
peratures have attracted much interest. To discuss giant 
resonances in hot nuclei, RPA has been extended by in- 
cluding thermal effects [10-17] . Such extension is referred 
to as finite temperature RPA or thermal RPA. 

Recently, there have been renewed interests in atomic 
nuclei at finite temperatures, in connection to nuclear 
astrophysics [18-22]. See also Ref. [23] for collective exci- 
tations in hot exotic nuclei. The finite temperature RPA 
has often been used to estimate e.g., beta decay rates in a 
stellar environment [19, 20, 24]. In order to make quan- 
titative calculations for nuclear astrophysical purposes, 
especially for r-process nucleosynthesis, it is necessary to 
know the accuracy of finite temperature RPA also at rel- 
atively low temperatures. 

The aim of this paper is to assess the applicability of 
finite temperature RPA using the Lipkin-Meshkov-Glick 
model [25]. This is a schematic solvable model and has 
been employed extensively to test many-body methods 
[1]. A similar study for finite temperature RPA has al- 
ready been done by Rossignoli and Ring [26] (see also Ref. 
[27]), but they explicitly investigated only the positive en- 
ergy part of a strength function. Here we investigate both 
the positive and negative energy parts separately. Since 
the sensitivity to the thermal occupation probability is 
large for deexcitation process at low temperatures, such 
investigation provides an interesting test of finite temper- 
ature RPA. We also consider both the canonical and the 
grand canonical ensembles for the exact solutions of the 
Lipkin model, while Ref. [26] considered only the grand 
canonical ensemble. This is important because the finite 
temperature RPA is based on the grand canonical ensem- 
ble despite that the number of particle is well-defined in 
actual nuclei. For compound nuclei formed in heavy- 
ion fusion reactions, the grand canonical ensemble may 
be justified because of neutron evaporation processes[14]. 
However, it is not obvious whether the same argument 



holds for nuclei in a stellar condition at low tempera- 
tures. By comparing the results with canonical ensem- 
ble to those with grand canonical ensemble, one can get 
some insight about the applicability of many-body theo- 
ries based on the grand canonical ensemble, such as finite 
temperature RPA. 

In the Lipkin model, one considers two single-particle 
levels, at energy of — e/2 and e/2, respectively, each of 
which has 2il-fold degeneracy. The Hamiltonian for this 
model reads [25] 

H = ek -^(K+K++,k-K-), (1) 

where V is the strength of a two-body interaction. The 
operators Kq, K+ and are defined as, 
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Here, and c]^ are the creation operators for the lower 
and upper levels, respectively. 

The exact solutions of the Lipkin model can be ob- 
tained with the quasi-spin formalism [1, 25]. The eigen 
states are then classified in terms of the eigen value of 

the operator K 2 = Kq + (K + K + + K-K-)/2. Denoting 
those states and their energy as | J a) and Ej a , respec- 
tively, the strength function for the canonical ensemble 
is given by [26] , 
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where /3 = 1/ kT is the inverse temperature and F is the 
transition operator. We have assumed that F does not 
change the value of J. Yc(J) is the degeneracy of the J 
state given by [28] 
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where N is the number of particle in the system, and 
J max = min[iV/2, 217 — N/2] is the maximum J for given 
N. Zc in Eq. (4) is the partition function given by 
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In this paper, we consider only a system with N = 20, 
(that is, half- filling) . In this case, the chemical poten- 
tial \x in the grand canonical ensemble is zero, and the 
exact strength function for the grand canonical ensemble 
is given by a similar formula as in Eq. (4) but with a 
different value of degeneracy [28], 



Y GC (J) = W GC (J) - W GC (J + 1)(1 - <fj, jM J, (8) 
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In addition to the exact solution, we also seek an ap- 
proximate solution for the strength function using the 
finite temperature RPA. To this end, we first solve the 
thermal Hartree-Fock equation[l, 10, 27] 
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is the thermal occupation probability of the Hartrcc-Fock 
state k. With the Hartree-Fock basis, we assume that the 
excitation operator of the system is given by 
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where and a^ oi are the creation operators for the 
Hartree-Fock states. The finite temperature RPA equa- 
tion then reads [10], 
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with 



a = e 1 -e + (f -f 1 )--(N-l)Vsin 2 2a, (17) 
b = -(/o - /i) • (N - 1)V sin 2 2a) , (18) 



where Doo = cos a. With the solutions of the RPA 
equations, the RPA strength function is obtained as 
[10, 12, 13, 29], 
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where the matrix elements for F are taken with the 
Hartree-Fock basis. 

Let us now solve the model Hamiltonian numerically 
and compute the strength function. For this purpose, 
we take the particle number to be A' = 20, and set 
VN/e — 0.5. As a transition operator, we consider 
F = (K+ + K-)/2. Following Ref. [26], we smear the 
strength function with a width of rj/e — 0.1. 

Figure 1(a) shows the strength function at tempera- 
ture of T/e=0.05. For our choice of parameters, the dif- 
ference in the strength function between the canonical 
and the grand canonical ensembles is small, and we only 
plot the result of grand canonical ensemble as the exact 
solution (see the solid line). The dashed and the dotted 
lines are the result of finite temperature RPA and ther- 
mal Hartree-Fock, respectively. At this low temperature, 
the thermal effect is almost negligible, and the strength 
function is actually almost the same as that at zero tem- 
perature. Notice that the RPA well reproduces the exact 
strength function. The Hartree-Fock result is not sat- 
isfactory, and the RPA correlation plays an important 
role. 

At finite temperatures, the excited levels are thermally 
occupied with a finite probability. The probability for de- 
excitation of the excited states with the operator F ap- 
pears in the negative energy part of the strength function. 
Fig. 2 shows the negative energy part of Fig. 1 (a). Since 
the temperature is low, the thermal occupation probabil- 
ity of the excited states is negligibly small. Nevertheless, 
we find that the RPA works reasonably well. The energy 
of the first excited state is 0.897e for the exact solution, 
while it is 0.88e in RPA. As the energy is slightly under- 
estimated in RPA, the peak of the strength function is 
somewhat overestimated. Despite this, we will show later 
that the temperature dependence of total strength is well 
reproduced with finite temperature RPA (see Figs. 3 and 

The strength functions at higher temperatures, 
T/e=0.5 and 1.0, are shown in Figs. 1(b) and 1(c). At 
these temperatures also, one sees that the finite tem- 
perature RPA works well both for the positive and the 
negative energy parts. Especially, the shift of the peak 
position in the strength function due to the finite temper- 
ature effects is well reproduced with RPA. As the temper- 
ature increases, even the thermal Hartree-Fock method 
reproduces the exact strength function. 

The total strength defined as 
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FIG. 1: (Color online) The strength function for the oper- 
ator F = (K+ + K-)/2 obtained with the several methods 
for VN/e=0.5. Figs. 1(a), 1(b), and 1(c) correspond to the 
temperature of T/e — 0.05, 0.5, and 1.0, respectively. The 
solid line is the exact result for the grand canonical ensemble, 
while the dashed and the dotted lines denote the solutions 
of finite temperature RPA and thermal Hartree-Fock, respec- 
tively. The exact result for the canonical ensemble is almost 
the same as the solid line, and is not shown in the figure. 



is plotted in the upper panel of Fig. 3 as a function of 
temperature. The lower panel shows the ratio of the total 
strength to that at zero temperature. The result of finite 
temperature RPA closely follows the exact result of grand 
canonical ensemble, as has been noted in Ref. [26]. We 
have confirmed that this conclusion remains qualitatively 
the same even for a smaller number of particle number, 
e.g., N = 10. The result of canonical ensemble, shown by 
the thin solid line, is close to the result of grand canonical 
ensemble, although the difference is not negligible. The 
thermal Hartree-Fock, on the other hand, leads to an 
inconsistent temperature dependence of total strength, 
as can be seen in the lower panel. 

The contribution of the negative energy part to the 
total strength is shown separately in Fig. 4. As one 
can see, the finite temperature RPA yields the correct 
temperature dependence of total strength even at very 
low temperatures. 
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FIG. 2: (Color online) Same as Fig. 1(a), but only with the 
negative energy part. 
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FIG. 3: (Color online) The total strength obtained with 
the several methods as a function of temperature (the up- 
per panel). The meaning of each line is the same as in Fig. 1, 
except for the thin solid line which denotes the exact result 
of canonical ensemble. The lower panel shows the ratio of the 
total strength to that at zero temperature. 



Let us next discuss briefly a case with a stronger 
coupling. Figure 5 shows the strength function for 
VN/e=2.0 at temperature T/e=0.675. This corresponds 
to Fig. 2(b) in Ref. [26]. In this case, the canonical 
and the grand canonical ensembles yield slightly differ- 
ent strength distributions from each other. That is, the 
grand canonical ensemble leads to a smoother strength 
function because of number fluctuation, although the 
overall behavior is similar to each other. One also notices 
that the finite temperature RPA significantly underesti- 
mates the thermal broadening of the strength function 
a/o the ground state correlation. However, as has been 
argued in Ref. [26] , the finite temperature RPA provides 
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FIG. 4: (Color online) Contribution of the negative energy 
part to the total strength shown in Fig. 3. The thin and the 
thick solid lines are indistinguishable on this scale. 



and the finite temperature RPA is satisfactory. We have 
checked that this is the case even for a stronger coupling, 
VN/e=4.0. 



TABLE I: The total strength for VN/e=2.Q at temperature 
T/e=0.675 obtained with several methods, that is, the ex- 
act result with canonical (C) and grand canonical (GC) en- 
sembles, thermal RPA (TRPA), and thermal Hartree-Fock 
(THF). The contribution from the negative energy part is 
also listed in the parentheses. 
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FIG. 5: (Color online) Same as Fig.l, but for a stronger cou- 
pling strength, VN/e=2.0 at temperature T/e=0.675. 



a reasonable estimate for the total strength. Table I sum- 
marizes the total strength for this particular choice for 
the parameters. The agreement between the exact results 



In summary, we have investigated the applicability 
of finite temperature RPA using a schematic solvable 
model. We have shown that the finite temperature RPA 
provides a reasonable estimate for the total strength, 
both for the excitation and the decay processes. This is 
the case even at low temperatures. For a small coupling 
case, the finite temperature RPA also yields a reason- 
able strength function itself. We have also shown that 
the canonical and the grand canonical ensembles lead to 
similar strength functions, as well as the total strengths, 
to each other. We thus conclude that the finite tempera- 
ture RPA, being based on the grand canonical ensemble, 
provides a reasonable tool to discuss properties of hot and 
warm nuclei, including those in a stellar environment. 
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